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Abstract: 

Dynamical systems exhibiting both PT and Supersymmetry are analyzed in a general 
scenario. It is found that, in an appropriate parameter domain, the ground state may or may not 
respect PT-symmetry. Interestingly, in the domain where PT-symmetry is not respected, two 
superpotentials give rise to one potential; whereas when the ground state respects PT, this 
correspondence is unique. In both scenarios, supersymmetry and shape-invariance are intact, 
through which one can obtain eigenfunctions and eigenstates exactly. Our procedure enables one 
to generate a host of complex potentials which are not PT-symmetric, and can be exactly solved. 

Introduction 

Complex potentials symmetric under simultaneous Parity (P) and Time -reversal (T) 
transformations are known to yield real energy eigenvalues in suitable parameter domain, 
sharing common eigenfunctions with the PT operator [1-3]. Beyond this range, these PT- 
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symmetric potentials correspond to complex-conjugate spectra with eigenfunctions connected 
through PT operation. This spontaneous breaking of PT-symmetry is expected as PT is a non- 
linear operator. Various models have been studied, both numerically and analytically, to 
illustrate the above structure [4, 5]. When these systems are analyzed under Supersymmetric 
quantum mechanics (SUSY-QM) [6], in the real domain, they are found to be isospectral to a 
real potential [7]. We here put forward an approach to arrive at both real and complex-conjugate 
domains starting from a class of superpotentials, with complex parameters. For the real case, it is 
possible to have a unique superpotential. For the complex-conjugate domain we have multiple 
superpotentials, each representing a separate section of the complete Hilbert space, but 
possessing the same potential. Moreover, these two classes of superpotentials are related through 
simple parametric variation, showing a bifurcation in the complex energy plane. 

In the following, we will discuss SUSY-QM briefly and then directly proceed to 
discuss a particular example of PT-symmetric potential [8]. Analysis of other complex potentials 
which are not PT-symmetric will be discussed, followed by some recent experimental findings. 

1. SUSY-QM: A Brief Introduction 

Following the factorization approach to analytically solvable Hamiltonians [9], in 
SUSY-QM [6, 10, 1 1, 12], the Hamiltonian can be written as H_(x) = A* A, where 4 + = 
— d/dx + W(x) and A = d/dx + W(x); with the superpotential W(x) being defined as: 
V + (x) = VK 2 (x) ± dW I dx. Here V_ is the potential corresponding to H_(x) and V + corresponds 
to another Hamiltonian H + (x) = AA^ , which is the superpartner of H_(x). The eigenstates of 
H ± (x) have one-to-one correspondence: 



0+(x) = [E- +1 ]- 1/2 Axp- +1 (x), 

and, 0" +1 (x) = [E +r 1,2 AHi(xy, (1.1) 

except for the ground state of H_(x), which is defined as Axp (x) = 0, and is expressible as 

xPodx) = e -S Xw ^ dx '. (1.2) 
Also, from (1.1), the energies are related as E£ — E~ +1 , whence the potentials V+(x) are called 
isospectral. Further, Gendenshtein [13] showed that, if two isospectral potentials are related as 
V+(x; a ) = V-(x; a x ) + R(a{), where a is a parameter in V+, % = /(a ) and i?(a x ) is 
independent of x; then the potentials are called shape-invariant. From such potentials, one can 
construct a hierarchy of Hamiltonians: 

H s = -d 2 /dx 2 +V_(x; a s )+Z s k=1 R(a k ) = - d 2 /dx 2 +V + (x; a s . t ) + Z|=i^(a fc ), 
with ground state energy = Zk=i R(. a k)- O n identifying, Z/ 1 = H + and H° = H_, the energy 
of the nth level of H_ is found to be Eq = Y,k=i ^( a /c)- Further, the excited states of the original 
Hamiltonian can be found from the ground state of H n as: 

0~(x; a s ) oc A^{x; a M + (x; a x ) i4 t (x; On-Jfo (x; a n ). (1.3) 

Therefore if the potential is shape-invariant, the eigenvalue problem can be solved 
completely by this algebraic method and the spectrum with corresponding eigenfunctions can be 
obtained. The application of the SUSY-QM approach to PT-symmetric pseudo-Hermitian 
Hamiltonians has been carried out in detail [14, 15], including the construction of the appropriate 
norm. Recently [16], the procedure to construct the norm for general pseudo-Hermitian 
Hamiltonians has been shown. 



2. Construction of Complex PT-symmetric Potentials 



Non-Hermitian, complex PT-symmetric potentials are known to have real eigenvalues over 
a definite range of parameters appearing in them, beyond which they show complex-conjugate 
(CC) spectra, with corresponding wavefunctions being PT-related. If such potentials are shape- 
invariant then the eigenvalue problem can be solved through SUSY-QM. 

We start with general superpotentials with complex parameters, leading to both real and 
CC spectra for different range of the parameter values. In the real domain, for a given potential, 
the superpotential is unique, for a given parameter range. For the CC domain, two different 
superpotentials yield the same potential, for a different range of parameterization. Thus, we 
arrive at the SUSY condition for phase-transition of the spectrum from real to CC values owing 
to spontaneous PT-symmetry breaking. 

For example, we will consider generalizations of PT-symmetric shape-invariant 
complex potentials, V(x; a) = — y 1 sec/i 2 ((2x) — iV 2 sech(ax) tanh(ax), which was analytically 
solved by Z. Ahmed [8]. Here V 12 and a are real constant parameters. We propose a 
superpotential, 

Wf T = {A± iC PT ) tanh(ax) + {±C PT + iB) sech(ax), (2.1) 
with A, B, a are real parameters. The corresponding potentials are, 

V±(x) = -[04 ± iC PT XA ± iC PT + a) - {+C PT + iB) 2 ]sech 2 {ax) 

-i{±C PT - B)[2{A ± iC PT ) + a] sech(ax) tanh(ax). (2.2) 
This potential is not PT-symmetric, for which the co-efficient of the even-P term has to be real 
and that of the odd-P part has to be imaginary. Invoking this condition, we arrive at the 
bifurcation condition, 

C PT [2{A -B) + a] = 0. (2.3) 



From (2.3), there can be two possibilities. For C = 0, we have from (2.1), 

W PT (x) = W real (x) = i4tanh(o:x) + iB sech(ax), (2.4) 

giving, from (2.2), 

V±(x) = V-(x) = —[A(A + a)+ B 2 ]sech 2 (ax) + iB(2A + a)sech (ax)tanh (ax). (2.5) 
Then, from (2.4), we get the real spectrum of the above potential through shape-invariance as, 

E = -(na-A) 2 , (2.6) 
modulo a constant term, and from (1.3), the eigenfunctions as [9], 

A A B 1 A B 1 
A \ B 1 > 1 

i/j n (x) cc [sec (ax)]a exp\—i-tan~ 1 {sinh(ax)}\P n a a 2 ' a a 2 [i sinh (ax)], (2.7) 
Again, if C PT 0, from (2.3), one gets A = B — ^, which when substituted in (2.1) and (2.2), 
yields, 

W,r T {x) = VK c ±(x) = 04 ± iC PT ) tanh(ax) + [±C PT + i(A + |)] sech (ax), (2.8) 
and, K_±(x) = K_ c (x) = - [2,404 + a) - 2{C PT ) 2 + ^] sech 2 {ax) + 

i [2404 + a) + 2{C PT ) 2 + sech(ax) tanh(ax). (2.9) 

Corresponding to two different superpotentials, we arrive at the CC spectrum, 

£± = 2n{A ± iC PT )a + (na) 2 , (2.10) 
modulo a constant term again, and the eigenfunctions as, 

oc [sech (ax)]^ A±lcPT) exp [{ - \{ A + ^)+ tan _1 {sin/i(ax)}] x 

P n ±t2 ~' 2 " + ^[isinh(ax)]. (2.11) 
The specific parameterization condition C PT ^ is the SUSY criterion for spontaneously broken 
PT, which is different from the analytic parameter criterion with PT being unbroken over a range 
of parameters. Further, for broken PT, two superpotentials corresponding to the 



bifurcation in the Hilbert space results into the same complex PT-symmetric potential. This non- 
uniqueness is different from isospectral deformation as it arises parametrically. In the complex 
case shape-invariance leads to complex energy shifts, while in the real case this shift is real, 
which was observed earlier [17]. The key structure of such potentials is that the odd Parity part 
has purely imaginary co-efficient and the even part has purely real co-efficient. More examples 
are listed in the source paper. 

Conclusion 

One can construct non-PT-symmetric shape-invariant complex potential through minimal 
complexification of the superpotential corresponding to the real counterparts []. For example, the 
Poschl-Teller potential [18] U(x) = U a sech 2 (ax) + U b csch 2 (ax), where U ab and a are 
constant parameters, can be complexified by considering the superpotentials 
W x — Atanh(ax) + iBcoth(ax) or W 2 = iAtanh(ax) + Bcoth(ax), A, B and a being 
constant parameters. In both cases, the spectra are complex. The non-uniqueness of the real 
spectra for the real potential is lifted by the complexification, with the imaginary part being 
equispaced. Further, the wave-functions become normalizable over a greater parameter range. A 

tec 

complex radial Coulomb potential can be constructed from the superpotential, W(r) = — + B, a 

and P being independent of r, which results into complex principal quantum numbers, with the 
wave-function being normalizable over a greater parameter range than those corresponding to the 
real counterpart [19]. Recently, spontaneous PT breaking has been observed experimentally in 
optical fibers, where the gain-loss profile showed bifurcation in the complex plane beyond a 



critical value of the optical loss co-efficient [20]. We propose similar observation in terms of our 
parameters, if they are identified properly for an optical system. 

# Recently it has been shown for unbroken PT-symmetry, the potential in Eq.(2.5) corresponds 
to not one, but two superpotentials, under the additional sl(2) symmetry of the system [21,22]. 
But each of them is shown to be independently mapped to the same pair of superpotentials in 
the broken PT sector under the SUSY parameterization [23]. 
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